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, 1 INTRODUCTION 



ABSTRACT 

Hardness-intensity diagrams (HIDs) have been used with great success to study the 
accretion states and their connection to radio jets in X-ray binaries (XRBs). The anal- 
ogy between XRBs and active galactic nuclei (AGN) suggests that similar diagrams 
may help to understand and identify accretion states in AGN and their connection 
to radio loudness. We construct "disc-fraction luminosity diagrams" (DFLDs) as a 
generalization of HIDs, which plot the intensity against the fraction of the disc contri- 
bution in the overall spectral energy distribution (SED). Using a sample of 4963 Sloan 
Digital Sky Survey (SDSS) quasars with ROSAT matches, we show empirically that 
an AGN is more likely to have a high radio:optical flux ratio when it has a high total 
luminosity or a large non-thermal contribution to the SED. We find that one has to 
consider at least two-dimensional diagrams to understand the radio loudness of AGN. 
To extend our DFLD to lower luminosities we also include a sample of low-luminosity 
AGN. Using a simulated population of XRBs we show that stellar and supermassive 
BHs populate similar regions in the DFLD and show similar radio/jet properties. This 
supports the idea the AGN and XRBs have the same accretion states and associated 
jet properties. 

Key words: accretion, accretion discs - black hole physics - galaxies: active - 
quasars:general - ISM: jets and outflows - X-rays:binaries 
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In spite of a vast difference of scales, both black-hole 
X-ray binaries (XRBs) and active galactic nuclei (AGN) 
arc thought to be powered by an essentially scale- 
invariant central engine consisting of a black hole, an 
accretion dis c surrounded by a corona and a relativis- 
tic jet fe.g.. Ish akura fc Sunvaevl Il973l : lAntonuccl Il993t 
iMirabel fe Rodriguez! Il999lh In the study of AGN, it is 
a long-standing puzzle that some objects show strong ra- 
dio emission while in others only little radio emission 
is found, even though their optical propert ies are simi- 
lar (iKellermann et al.lll989t iMiller et al.lll993h . XRBs, i.e., 
stellar-mass accreting black holes, show states associated 
with strong radio emission as well as states during which 



jTananbaum et al. 
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). The 



have been ascribed t o different jet prope rties of the different 
accretion states fe.g.. lFender et all2004!) . Due to the shorter 
timescales of accretion onto XRBs, it is possible to study the 
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different accretion states and their association with jets in 
detail. 

The quantitative comparison of XRBs and AGN 
has made a major step forwards in recent years with 
the dis covery of a fundamental plane of blac k-hole ac- 
tivity iMerloni et all 120031: iFalcke et all 12004!) . This is 
a correlation between radio and X-ray luminosities 
and black hole mass which fits both classes of ob- 
jects (both individually and as an ensemble). Also the 
variability p roperties of AGN and XRBs show many 
similarities JUttlev et al] 120021: iMarkowitz et all 12003: 
iKording fc Falckd 12004 lAbramowicz et al.l 12004) . There- 
fore, it is a promising route to use the concepts and methods 
developed to study the evolution of radio jets in XRBs for 
exploring the puzzle of the radio loudness in AGN. 

XRB outbursts seem to occur in cycles of several dis- 
tinct states that are si milar in every l arge outburst and fr om 
object to object (e.g., lNowaklll995l : iBelloni et alj l2005l see 
Fig. 0. At the beginning of the cycle the source is in the 
hard state characterised by a hard power law in the X- 
ray spectrum . At radio fr equencies, one ususally observes 
a steady jet (Fendei||2 p^)l|h In this sta te the standard ac- 
cretion disc 1 Shakura fc Sunvaevl Il973h is probably trun- 
cated or unobservable for some other reason and the inner 
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part of the accretion disc replaced by an inefficient accretion 
flow (e.g-, an advection dominated accretion flow s, AD AF: 

aT1 ll997l but 



iNaravan fc Villi 9941: 


Chen et a.Ull995l 


see also iMiller et al. 


2006). The X-r 



ll975l : ISunvaev fc Trunrperll97flf ). However, it is also possible 
that the compact jet may contribute to the X-r ay emission 
or even dominate it (iMarkoff et all 1200 ll Eooj) . Once the 
source leaves the hard state, it is often found in an inter- 
mediate state (IMS). One finds two IMSs: first, the source 
enters a hard IMS characterised by a hard spectral compo- 
nent and band-limited noise in the power spectrum. Also in 
this state a jet is observed in the radio. Later, a soft IMS is 
found, dominated by a soft spectral co mponent with power - 
law noise in the power spectrum (e.g.. Eelloni et aT| [2005'). 
Approximately at the transition from the soft to the hard 
IMS one usually finds a radio flare and the ejecta can of- 
ten be resolved in the radio band. In the soft IMS no radio 
emission is found. After leaving the IMS, the source may go 
to the soft state, where the X-ray spectrum is dominated by 
a soft thermal component that is thought to originate from 
a standard accretion disc. The radio jet is still q uenched in 
this state jFender et al.lll999l : ICorbel et al]|2000T . 

It has not been established definitely how the phe- 
nomenology of AGN can be mapped onto these XRB states. 
Quasars are AGN with bright optical continuum emission, 
which are thought t o have high accretion r ates and strong 
disc emission (e.g., lUrrv fc Padovanl Il995l) . Here, we use 
the term AGN in a general sense to mean all supermas- 
sive accreting black holes, and we use the term "quasar" 
irrespective of the radio properties, i.e., synonymous with 
"quasi-stellar object" (QSO), and without imposing any lu- 
minosity cut. 

For AGN, one commonly adopted definition of a "radio- 
loud" source is for the R parameter, the ratio of the radio 
flux (at 1.4 GHz) to the optical B-band flux, to be larger 
than 10. We follow this convention here. It has been known 
for over two decades that radio-loud quas ars have h igher 



X-ray fluxes that radio-quiet objects (e.g., Zamorani et alJ 
ll98ll:lElvis et alJll994l : IZdziarski et al.lll995l) . perhaps hint- 



ing (with hindsight) that they represent different accre- 
tion states. iPounds et al.l lll995l) suggested that narrow-line 
Seyfert 1 objects may be analogues of the soft state. The 
fundamental plane of accreting black holes has been inter- 
preted as evidence that low-luminosity AGN (LLAGN) ar e 
hard-state objects llFalcke et alJl2q04l:lK5rding et al]l200rj) . 
This has also been suggested bv lHdl2005li ~"based on the lack 
of an observable "big blue bump" in the spectral energy dis- 
tributions (SE Ds) of LLAGN. Again using the fundamental 
plane. iMaccarone et alJ i2003l) suggested that AGN with an 
intermediate accretion rate (a few percent of the Eddington 
rate) are the analogue o f soft-state XRB s. 

On the other hand, iBorosorJ (120021) has suggested that 
radio-loud and radio-quiet quasars occupy distinct regions 
of a two-dimensional diagram of physical acc retion rat e M 
against Eddington-scaled luminosity L I L-pm . iBorosorj con- 
siders these two quantities a s the likely physical ori gin of 
the "eigenvectors" 1 and 2 (iBoroson fc Greer] I1992T) . and 
eigenvector 1 appears to be related to the radio proper- 
ties. The divisio n between radio -loud and radio-quiet ob- 
jects proposed bv lBorosonl (120021) i s one in black-hole mass, 
similar to those reported elsewhere iFranceschini et alll998l : 



[ Lacy et al.ll200ll: Ijarvis fc McLurell2002l) . However, neither 
tHTj^OO^T nor I Woo^^Uravri200^r found any evidence for a 
mass threshold between radio-loud and radio-quiet objects. 
Thus, the physical mechanism making some sources much 
more radio-loud than others remains a puzzle. 

For X-ray binaries, hardness-intensity diagrams (HIDs) 
have been used with great success to study the evolution 
of outbursts and separate the different accretion states (e. ; 
iHoman et allbOOlt iFender et all 12004 IBelloni et all 12005 
The radio-loud and radio-quiet phases of a source are well 
separated on this HID. Additionally, the timing properties 
of the X-ray lightcu rve change dramatically with the posi- 
tion in the HID (e.g.. lHoman et al-^OOil IBelloni et al.l2005l : 
lRemilrardl F2005). Hence, the construction of equivalent dia- 
grams promises to be a useful tool for the study of AGN, too. 
For XRBs, X-ray HIDs are useful diagnostic tools because 
spectral features of both the accretion disc and the hard 
power-law component are visible in the X-rays. The black- 
body temperature of the accretion disc scales with black 
hole mass as M _1//4 so that the SEDs of accretion discs 
around supermassive black holes peak in the optical/ultra- 
violet and an X-ray HID for AGN would not contain any 
information about the disc emission. Therefore, we have to 
construct a new type of diagram which is applicable to the 
study of XRBs as well as of AGN. In this paper, we present 
the definition of an analogous diagram for AGN (S|5J and use 
our method to explore the origin of radio loudness in a sam- 
ple of quasars from the Sloan Digital Sky Survey (SDSS) 
and LLAGN (©, also considering the impact of selection 
effects and other biases O- We compare the behaviour of 
AGN and X-ray binaries in JSJand conclude in SJH] 



2 METHOD 

2.1 Scaling hardness flux diagrams from XRBs to 
AGN 

X-ray HIDs are used to study the evolution of an XRB dur- 
ing an outburst dFender et alJl2004l : IBelloni et alJl2005l) . On 
the ordinate (y-axis) of these diagrams, one usually plots 
the X-ray count rate, while the abscissa (x-axis) denotes 
the hardness ratio (e.g., 6.3-10.5 keV count rate/3.8-6.3 keV 
count rate). For a sketch of a typical HID see Fig.Q While a 
source is in the hard state, its hard power-law, generally as- 
sumed to arise via Comptonization, dominates the spectrum 
and the source is on the right-hand (hard) side of the dia- 
gram. The X-ray emission in the soft thermally- dominated 
state is believed to be mainly due to a standard thin accre- 
tion disc (multi-temperature black-body). In this state, the 
source moves towards the left of the diagram. In AGN, the 
thermal emission of the accretion disc peaks in the optical 
or ultra-violet — far from the observed X-ray bands. Thus, 
an HID based solely on the X-ray spectrum cannot probe 
whether an AGN is in a "soft" or "hard" state. 

An HID that can be used to diagnose AGN as well as 
XRBs has to be based on parameters that are independent of 
the observing band, as at least the temperature of the accre- 
tion disc and therefore the peak of the emission component 
scales with black hole mass. The location of an individual 
source in the HID is determined by the total luminosity 
of the system and the relative strength of the power law 
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Figure 1. Sketch of an HID found for XRBs. For the abscissa 
one usually uses the hardness ratio and the ordinate represents 
the measured X-ray counts, i.e., the total measured X-ray flux. 
The arrows indicate the movement of a source during outburst. 



component compared to the disc component of a source. In 
XRBs, both the disc and the power-law component are ob- 
servable in the X-ray band. The X-ray luminosity is (caveat 
bolometric corrections) equal to the luminosity Lu + Lpl, 
where Lu denotes the disc luminosity and Lpl the lumi- 
nosity of the power-law component. The hardness reflects 
which component dominates the spectrum. Thus, jT^Z~ 
has a similar behaviour as the hardness ratio. We will refer 
to this quantity as the non-thermal fraction. This value ap- 
proaches zero for disc-dominated sources and 1 if the power 
law is the dominating component. 

In AGN, the disc luminosity can be measured in the 
optical band for non-obscured sources, while X-ray observa- 
tions provide information only about the power-law compo- 
nent. A benefit of the term , h f\ — is that it stays finite 
if one of the components approaches zero, which happens 
when a source is in the hard or the soft state. As it is a 
ratio of luminosities, it is furthermore independent of the 
distance. For the case that the luminosities of both spectral 
components scale similarly with BH mass, it is also inde- 
pendent of the BH mass. We will refer to diagrams plotting 
Lu + ipL against L -j^ — as disc-fraction/luminosity dia- 
grams (DFLDs). 

We note in passing that the ratio L L ^ — is proba- 
bly connected to the properties of the accretion disk wind 
at least in AGN. The UV luminosity (Lu) is a measure of 
the radiation pressure while the X-ray luminosity (Lpl ) in- 
fluences the ionisation state (see, e.g., iRicharda I2006I . and 
references therein). 

For XRBs we can observe the evolution of an entire 
outburst of an individual source and obtain a closed loop on 
the HID or DFLD, as the t imescales of an outburst is months 
to years llChen et alll997l) . For each position in the diagram 
we can observe whether or not the jet is visible. However, 
for AGN it is impossible to observe a full outburst cycle 
(quiescent — > quasar phase — + quiescent) as it would last for 



hundreds of millions of years (assuming a linear scaling of 
timescales with BH mass). Instead, we have to study a large 
population of AGN in the DFLD and observe their average 
radio loudness (i.e., the radio:optical flux ratio) as a tracer 
of the jet activity. 



2.2 Our sample 

To construct a DFLD for AGN we need a sample of AGN 
with well measured optical (Lu), X-ray (£pl) and radio 
fluxes. 



2.2.1 SDSS Quasars 

As our main sample, we begin with the lists of all objects 
from the spectroscopic dat abase 1 of SDSS DR5 (J. Adelman - 
McCarthy et al., in prep.; lAdelman-McCarthv et"afll200dl 
that have been identified as a quasar with high confidence. 
The sample includes all objects from the statistically com- 
plete quasar survey feichards et aljf2002h as well as quasar 
spectra selected in other ways, e.g., candidate optical coun - 
terparts of known X-ray sources jAnderson et al.ll2003h . 
Given the criteria used by the SDSS to classify objects as 
quasars, all objects are Type 1 quasars with strong broad 
emission lines (and potentially broad absorption lines). Ap- 
pendix EI gives the SQL query necessary to reproduce our 
initial sample. For X-ray and radio measurements, we use 
data from th e RQSAT All-Sky Survey (RASS, energy band 
0.1-2.4 keV; IVoges et al"lll999fl and the Very Large Array 
(VLA) Faint Images of the Radio Sky at Twenty centime- 
ters (FIRST) survey JWhite et alJl997l observing frequency 

1.4 GHz;) as given in the SDSS database. The radio fluxes 
are intended to be a measure of the core radio flux of our 
sources (see H4.41 . We only include quasars in the redshift 
range 0.2 ^ z ^ 2.5. Below a redshift of 0.2, not many 
quasars are found, and spectra of low-luminosity AGN are 
often dominated by galaxy light, making them more likely to 
be classified as galaxy by the SDSS pipeline. Redshifts above 

2.5 contribute only a small fraction of the total number of 
SDSS quasars, and not many of them have a ROSAT detec- 
tion at such high redshifts. In the adopted redshift range, 
there are a total of 64248 objects of which 5174 have a ra- 
dio detection, 4963 have an X-ray detection, and 730 have 
detections in both bands. 

We apply a ./^-correction to the flux measurements in 
all three wavebands in order to remove apparent luminos- 
ity differences that are caused by observing objects at dif- 
ferent redshifts, and therefore at different rest wavelengths. 
The Jf-correction is conventionally de fined as mjntrinsic = 
^.observed - K(z) (see lHogg et~alll200l e.g.). For the SDSS 
quasars, we use reddening-corrected i-band PSF magni- 
tudes (identical to those used to flux-limit the SDSS main 
quasar sample) and a pply the i-band i^-correction given by 
iRichards etail (|2006L Table 4). These authors remove the 
average contribution of emission lines as function of redshift 
to obtain a pure continuum luminosity with an assumed 
power-law spectrum f v oc v~ ' 5 and define the correction to 
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Figure 2. Left: DFLD of 4963 SDSS quasars with ROSAT detections in the redshift range 0.2 $C z $C 2.5, showing the lower limit to the 
average radio loudness R as calculated from FIRST data. The contour levels are separated by 0.5 dex in R. Right: Upper limit to the 
average R. The contour levels are the same as for the lower limit. 



be at z = 2. The continuum /("-correction is easy to com- 
pute since the SDSS magnitude system is (nearly 2 ) an AB 
system (denned as mag = — 2.51og( . /W3631 Jy) for a sourc e 
with constant f. JOke fc Gunnll983tlAbazauan et 81.1 12003): 
the correction is given by K(z) — —2.5(1 + a) log(l + z) for 
/„ oc v a . We convert the if -correction to Mi(z = 2) back to 
the conventional Mi(z — 0) by adding the offset 0.596 ap- 
propriate for an object at redshift z — 2 wi th the assumed 
powe r- law spectrum with a = —0.5 (see iRichards et alJ 
2006, eqn. [1]). We then transform to the B-band using 
B — i — 0.3 as appropriate for the same power-law spectrum. 
Thus, the B-band magnitude is obtained from the SDSS i- 
band PSF magnitude, galactic extinction Ai as given in the 
SDSS database, and the combined emission-line and contin- 
uum -correction K opt as 



B = i + 0.3- Ai 



K, 



opt 



0.596. 



(1) 



However, we do not apply any correction for contributions 
from host galaxy starlight. This correction is most important 
for the lowest quasar luminosities and for the lowest-redshift 
objects. We discuss the impact of this simplification in H4.3I 

For the /("-correction of the radio fluxes we assume op- 
tically thin spectra with a spectral shape /„ cx u~°-\ i.e., 
we assume the same spectral shape for the radio and opti- 
cal continua. For the X-ray data, we assume a photon in- 
dex r = 2, i.e., a spectral shape f v oc v^ 1 (found to be 
appropriat e for both radio-loud and radio-quiet quasars by 
iGalbiati et afll200Sh which has /("-correction because the 
effects of bandwidth narrowing and observing a steep spec- 
trum at lower rest-frame frequencies cancel. 

As the SDSS quasar sample includes objects with a large 
range of black-hole masses, it is in principle necessary to ap- 
ply a mass correction to the observed luminosities. Accreting 
black holes that are observed as a quasa r all have si milar Ed- 
dington ratios, in the range 0.05-1 (e.g.. lJesterl2005l) . Hence, 
to zeroth order, the luminosity of a quasar is expected to 
scale as the Eddington luminosity, i.e., linearly with black 



hole mass. However, the black hole masses computed for 
12245 quasars from SDSS DR1 bv lMcLure fc Dunlod (120041) 
have a mean of (4.75 ± 0.44) x 10 s Mq, i.e., the intrinsic 
scatter in the black hole mass of quasars accessible to the 
SDSS spectroscopic survey is comparable to the scatter in 
the black-hole ma ss determination ( estimated to be as large 
as 0.5-0.6 dex bv IVestergaardir2004h . Hence, the black-hole 
masses of all SDSS quasars are consistent with being equal 
to the mean and we apply no mass scaling. This needs to be 
borne in mind when comparing this quasar sample to AGN 
samples with substantially different black hole masses. 

We use a cosmology with Qm = 0.3,0a = 0.7, Ho = 
71 kms^Mpc -1 . 



2.2.2 Low-luminosity sample 

While our main results will be based on the SDSS quasars, 
w e construct a lo w-luminosity sample based on the sample 
of lHo et alJ <ll997f) as a co mparison. We have obtained X -ray 
fluxes for this sample from lTerashima fc Wilso'r] <l2003l) and 
the followin g surveys in order of p reference: The Chandra 
v3 pipeline Iptak fc Griffiths! I2003P'). th e XMM serendipi- 
tous X-ray survey (Bar cons et al J 120021) . and the R OSAT 
HRI pointed catalog iROSAT Scientific Teaml |2000F). Ra- 
dio flu xes for the sample have been taken from Nagar et alJ 
(2005). For low- luminosity AGN, the optical emission is of- 
ten dominated by galaxy light. Therefore, we use a corre- 
lation between H/3 and the B-band magnitude llHo fc Pend 
l200ll): 



(2) 



logL H /3 = -0.34M S + 35.1 



to compute the optical absolute magnitude. This correlation 
seems to hold for quasars as well as Seyfert galaxies. We use 
the correlation also for our low-ionization nuclear emission 
region (LINER) objects even though the correlation has only 
been shown to hold for Seyfert objects and PG quasars. 



2 See the discussion at http : //www. sdss . org/dr5/algorithms/f luxcal ."iifaitlt p : //www. xassist . org 
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The average BH mass of the SDSS quasars is signif- 
icantly larger than that of the nearby LLAGN. Thus, for 
LLAGN it is important to consider their BH masses. The 
BH masses of the LLAGN sample are calculated from the 
M-cr relation iMerritt fc Ferra rcsc 2001 ) using velocity dis - 
persions from the Hypercat catalog {Prugniel et al.lll99£ft . 
We scale the luminosities linearly with mass to the average 
mass of the SDSS quasars. 



2.3 Calculating the component luminosities 

To compute the accretion-disc luminosity, we use the rela- 
tion between the absolute optical _B-b and magnitud e and t he 
approximate disc luminosity given bv lFalcke et"ai]lll995l) : 

log L D = -0AM B + 35.9 (3) 

This disc luminosity measure has a scatter of 0.29 dex. 

The power-law component Lpl is estimated from the 
X-ray luminosity. To obtain the power-law luminosity, we 
extrapolate the measured ROSAT counts to a flux in the 
0.5-10 keV band using a photon index of F = 2. For the hy- 
drogen column density we use an aver age Galactic column 
densi ty of uh ~ 1 X 10 21 cm -2 (e.g., iDickev fc Lockmanl 
1990). The hydrogen column density in the Quasar is usu- 
ally negligible of broad line Qu asars (for the low abs orption 
of broad line objects see, e.g., iRichards et~al]|2003h . addi- 
tionally, the large redshift reduces its effect. This yields a 
conversion factor of lcps = 3.1 x 10 -11 erg s _1 cm~ 2 . 
Where Chandra observations are available, we can measure 
the flux in this energy range directly. The luminosity of the 
power-law component is chosen to be 

Lpl = 3£ 5 _ loke y ( 4 ) 

This "bolometric correction" of a factor 3 is fairly arbitrary 
and was chosen to distribute the quasars uniformly over the 
a;-axis of the DFLD. For the purpose of our analysis, it does 
not have to correspond directly to the real bolometric cor- 
rection of the quasar, but is probably not far off. However, 
a different value creates only a slightly different mapping of 
the sources on the defined axis. This is similar to the defini- 
tion of the hardness ratio in the HIDs, which is also arbitrary 
and chosen to give a good distribution of the sources in the 
HID. Combining these "bolometric" and the ^-corrections 
from H2.2.1I into the factors c op t and cx, our non-thermal 
fraction is related to the observed fluxes in both wavebands 
by 1/(1 + c pt /opt /[cx/x]). 

To compare the radio loudness of the AGN as function 
of their position in the DFLD, we also calculate the param- 
eter R, which is defined as the ratio of the radio flux at 1.4 
GHz divided by optical B-band flux. For inclusion in this 
diagram, we only accept sources which have optical and X- 
ray fluxes, but we include sources with only an upper limit 
for R. Including the expressions for our if-correction terms 
from t|2.2.1l R is related to fluxes from the catalogues by the 
expression 

log .R = 0.4(5-*) (5) 

favezic et al]l2002ft . where t is the A'-corrected 1.4 GHz ra- 
dio AB magnitude given by t = — 2.5[log(/FmsT/3631 Jy) + 
0.51og(l + z)] and B is the isT-corrected £?-band magnitude 
from eqn. Q. 



3 RESULTS 

We can now use our samples of AGN to construct DFLDs, 
i.e., the {/-axis of this diagram denotes Ld+^pl, while on the 
a:-axis we plot Lj ^+l ■ F° r every source, we can calculate 
the position in the DFLD and its radio loudness R, or an 
upper limit to R. To study the dependence of radio loudness 
on location in the DFLD, we compute the average radio 
loudness of AGN at each position. We do so by binning 
the data in 10 x 10 cells of the DFLD and calculating the 
arithmetic mean of R for each cell. We only show those cells 
in the diagram that contain more than 10 sources. We treat 
the upper limits to R in two ways: First we set R to zero for 
all sources without a radio detection. This average i?i ow is a 
lower limit to the real average of R. Secondly, we set R to the 
upper limit of the corresponding radio measurement. This 
average, i?high, is an upper limit to the real average R. Thus, 
the real average R has to lie between these two values. We 
will use these values in our figures to detect selection effects 
caused by the FIRST survey's radio flux limit. 



3.1 SDSS Quasars 

In Fig. |5]we show both upper and lower limit of the average 
R for our 4963 SDSS quasars with ROSAT detections. The 
average R value is largest in the upper right corner and 
declines gradually towards the lower left corner. The average 
log R in the top right corner is 2.6 dex higher than in the 
bottom left corner, which corresponds to a factor 398 in R. 
The gradient is seen both in the upper and lower limit of R 
and the difference between both maps is very small. Thus, 
we conclude that this effect is also found in the real averaged 
R. The effect can be described as combination of two trends: 
First, the most radio-loud objects for a given disc luminosity 
(i.e., optical luminosity) are among those sources that are 
the most X-ray bright (see Fig. Additionally, brighter 
objects tend to be more radio-loud for a given disc fraction 
(a fixed a;- value in the DFLD). The first effect is already 
known as r adio-loud objects are known to be more X-ray 
bright fe.g.. lElvis et~ai]ll994tls"hen et al.ll200rj) . 

As the average has been computed separately for every 
cell, two different cells are statistically independent. Thus, 
even though the average in each cell does not have to be 
significant, the fact that the map is smooth and has a well- 
defined gradient in R indicates that this effect is not a sta- 
tistical artifact. To verify this observation mathematically, 
we have performed a Kolmogorov-Smirnov test of the dif- 
ferences between the distributions of radio loudness in two 
well-separated areas of the DFLD. We compare all sources 
within a "high-luminosity" area given by 0.3 ^ x ^ 0.5 and 
46 y ^ 47 to those in a "low-luminosity" area given by 
0.3 < x < 0.5 and 44 < y < 45.5. In Fig. we show the 
Log N - Log R diagrams for both areas. Already visual in- 
spection suggests that the "high-luminosity" area contains 
a significantly larger fraction of radio-loud objects. This is 
confirmed by the Kolmogorov-Smirnov test, which rules out 
that the two underlying probability distributions are identi- 
cal with a confidence level of 99.994 % (D = 0.19). We also 
verified that the radio-loudness distribution is different for 
areas with similar luminosity but different disc-to-power law 
ratios. 

To verify that the X-ray emission does have an effect on 
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Figure 5. Limits to average radio loudness and percentage of radio-loud sources (R > 10) as function of disc luminosity. Left: all 64248 
SDSS quasars in the rcdshift range 0.2 $C z $J 2.5. Right: only the 4963 X-ray detected sources. The error bars give the uncertainty of 
the mean in each of the bins. 
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Figure 3. Modified DFLD showing the average radio-loudness. 
The y-axis gives the optical disc luminosity Ljj instead of the 
sum of the disc and the power law luminosity in the previous 
maps. The most radio-loud sources are those with a large power 
law component in the spectrum. 



the radio loudness, we plot the average radio loudness as a 
function of the disc luminosity (i.e., the optical luminosity) 
in the left-hand panel of Fig. for all 64248 SDSS objects, 
independently of their X-ray detection. The average R has 
been estimated independently in 20 logarithmic luminosity 
bins. This plot of the radio loudness as a function of lumi- 
nosity can be compared to the right-hand panel of the same 
figure, where we show the dependence of the radio loudness 
on the disc luminosity for the 4963 X-ray detected sources 
only. While there is the clear trend visible in the X-ray de- 
tected quasars that brighter objects are more radio-loud, 
only marginal effects can be seen in the full sample. We 
show the upper and lower limit to the radio loudness as well 
as the fraction of radio-loud objects as defined by R > 10. 
The trend can be seen in all three quantities. We can also 
find a clear trend in the radio-loudness if we project the 



Figure 4. Log N-log R diagram for two rectangular areas in the 
DFLD diagram. The position and extent of the upper and lower 
box is given in the text. For both areas we first find all objects 
with a position on the DFLD in that area and construct log TV-log 
R diagrams for both areas individually. Solid lines indicate upper 
limits to the average _R while dashed lines represent its lower limit. 



DFLD onto an axis roughly corresponding to the diagonal. 
In Fig. © we show such a projection for the SDSS sample 
( compare IShen et al.ll2006l . Fig. 11). 



3.2 Hard-state objects: LLAGN 

So far, we have only considered the SDSS quasars. All 
quasars are thought to be strongly accreting objects with a 
standard geometrically thin, optically thick disc. To obtain 
information about the shape of the DFLD at lower luminosi- 
ties, which may be in a different, radiativ ely inefficient ac - 
cretion mode, we now include the sample of lHo et alJ lll997t) . 
Due to the small number of sources we now show every cell 
for this sample, even if it includes only one source. Hence, 
for the LLAGN, the plotted "mean" R may just be a single 
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Figure 6. Projection parallel to the black lines in the DFLD 
shown in the inset (with the chosen scaling of the coordinate axes, 
the line onto which objects are projected is not perpendicular to 
these lines). We project on to the line defined by A = (0.93 e x + 
0.36 e y ). This projection includes only sources with measured X- 
ray fluxes. 
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Figure 7. DFLD showing the average r adio loudness for SDSS 
quasars and LLAGN from the iHol ll999Fl sample. Note that the 
gap between LLAGN and Quasars is an artefact of our sample 
selection. 

measurement that could be substantially different from the 
true mean R. 

In Fig.Qwe show the joint DFLD for the SDSS quasars 
and the LLAGN sample. Not surprisingly, the LLAGN lie at 
lower luminosities than the SDSS quasars. However, except 
for a few outliers they all lie a t the right sight of the DFLD. 
This is also not surprising, as lHol l)l999f) already noted that 
LLAGN are less optically bright for a given X-ray luminosity 
than strongly accreting AGN (Seyfert galaxies and quasars). 
Furthermore, the LLAGN in our DFLD seem to be more 
radio-loud on average than SDSS quasars. As mentioned, we 
have to plot every cell for the LLAGN even if it only contains 
one object, so the mean value of R for an y given cell may not 
be close to the true mean. But already iHo fc Pend (1200 ll) 



have shown that the average low-luminosity Seyfert galaxy 
is more radio-loud than the average "radio-loud" PG quasar 
if only the nuclear luminosity is considered (i.e., removing 
the contribution of starlight). Thus, we find that LLAGN 
populate the bottom right-hand side of the DFLD and are 
— on average — more radio-loud than the SDSS quasars. 



4 CRITICAL ANALYSIS OF SELECTION 

EFFECTS AND ASTROPHYSICAL BIASES 

In this section, we consider whether the shape of the distri- 
bution of SDSS quasars in the DFLD and their average radio 
loudness could be caused entirely due to selection effects, or 
underlying correlations between the observables and other 
parameters. We consider in detail the selection effects of 
flux- limited samples ( H4.1H . black-hole mass scalings ( t|4.2H . 
orientation, beaming and host-galaxy light f <J4. 3B and the 
size- luminosity evolution of radio sources( i|4.4t . Based on 
our analysis, none of these effects can produce the entire 
signal we report. 

4.1 Selection effects 

Our main sample is basically a single quasar survey with 
a fairly bright flux limit — most quasars in our sample 
are drawn from that part of the quasar sample limited by 
i < 19.1. Even though we include several other SDSS tar- 
get categories, there is an implicit flux limit i < 20 imposed 
by the signal-to-noise that can be achieved in the exposure 
time of the spectroscopic observations in the SDSS. As the 
optical luminosity function of quasars is rather steep, this 
means that the dominant selection effect is that we only 
have a rather restricted dynamic range in quasar luminos- 
ity at fixed redshift, of order 10-30 (the SDSS spectroscopic 
survey has a bright limit i > 15 and substantial numbers 
of objects are only found at i > 16). Thus, only the most 
luminous object at any given redshift are found, and the 
variation in luminosity is largely due to a variation in red- 
shift. Thus, on average there is a clear trend that the redshift 
of the lowest luminosity sources is small (z — 0.2) while the 
highest luminosities are found at large distances (z = 2.5). 

To assess the impact of this well-known degeneracy be- 
tween redshift and luminosity in flux-limited surveys, we 
have verified that the effect is visible in small redshift slices, 
e.g., in the subsample with redshifts between z = 0.3 and 
z = 0.4 as shown in Fig. |H| Furthermore, we have verified 
that the ^-correction is not creating the observed trends in 
the DFLD. Even if we assume a flat radio spectrum or do 
not isT-correct at all, we find the observed effect. The differ- 
ence in radio loudness between the average SDSS quasar at 
low redshift (z ~ 0.2 ) and those at large redshifts (z fa 2.5) 
is 0.6 in log R. This difference is probably mainly due to a 
difference in accretion rates. However, even if this is purely 
due to evolution, we note that the observed range in R in 
the DFLD is around 2.5 dex, so evolution would only be a 
minor effect. This suggests that the radio loudness gradient 
in the DFLD is not an artifact of redshift evolution, but we 
cannot rule out that evolution plays a minor role. It clearly 
will be desirable to construct a DFLD with a sample that 
is not affected by evolution (i.e., one with a narrow redshift 
distribution) . 
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Figure 8. DFLD of the redshift slice 0.3 < z < 0.4. In this 
redshift range our sample contains 676 sources with optical and 
X-ray detections of which 70 also have radio detections. Due to 
the small number of available sources we show all bins with more 
than 5 sources. Similar to the DFLD of the full sample sources in 
the upper right corner have large average R values, while those 
in the lower left are less radio-loud. 



We have shown the average radio loudness as a func- 
tion of the position in the DFLD in two versions: as an 
upper limit and as a lower limit to the real radio loudness 
(i.e., the value we would observe with a radio telescope of 
infinite sensitivity) . More sensitive radio observations would 
increase the lower limits while they would decrease the up- 
per limit map. As there are only minor differences between 
both maps, we conclude that we our result is not affected 
strongly by the radio detection limit. 

As an additional test, we have performed a Monte- 
Carlo simulation to check whether the trends observed in 
the DFLD for the SDSS quasars are merely due to selection 
effects arising from the use of flux-limited surveys. As null 
hypothesis, we assume that radio, optical and X-ray fluxes 
of all objects are uncorrelated, so that the conditional lumi- 
nosity function is simply the product of the individual lumi- 
nosity functions: $(Lr,Lo,£x) = $i(£r) $a(£o) $3(Xx)- 
At the high luminosities that are accessible to our three 
surveys, all three luminosity functions are well-described by 
simple power laws oc j^R. o.x. For the radio luminos- 

ity function we use /3r = 0.78 jNaear et al.ll2005l) . for the 
X-rays fix = 2.34 llJeda et al.ll200J) and for the o ptical 
luminosity function f3o = 2.95 iRichards et alj|2006f) . For 
this simple test we assume that the power law indices do 
not change with redshift. The number of sources with radio, 
optical and X-ray detections at a given redshift is fixed to 
the observed value. This incorporates density or luminosity 
evolution with redshift, as well as the redshift- dependent se- 
lection efficiency of the SDSS quasar survey. The resulting 
DFLD of the simulated sample is shown in Fig. [5] The trend 
visible in the observed DFLD (Fig. [21 , i.e., an increase in the 
mean radio loudness from the bottom left to the top right, is 
not observed. Futhermore, we have repeated this luminosity 
function test with a population of artificially beamed steep- 
spectrum sources: a population of st eep-spectrum sources 
with (3b. = 2.31 iWillott et all 1200 ll) . 10% of which have 
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Figure 9. Monte-Carlo simulation of the radio loudness as a 
function of the position on the DFLD, assuming no correlation 
between radio, optical and X-ray luminosities. We assume simple 
power-law luminosity functions and give the simulated sample the 
same redshift distribution as our observed sample. 



their radio flux increased by a factor 100, but still without 
correlation to optical or X-ray properties. This simulation 
cannot reproduce our findings, either. 



4.2 Black-hole mass 

We are not able to correct the DFLD of our SDSS quasars 
for the BH mass as we do not know the mass for most of 
our objects. Even though most quasars in the SDS S have 
masses around 8.7 x 10 s M© (iMcLure fc Dunkvcfeooll . more 
luminous objects are more likely to have a more massive 
black hole. Thus, objects in the lower-luminosity part of 
the diagram will have a lower mass on average, while high- 
luminosity objects will have larger BHs. It has been sug- 
gested that the radio power of AGN correlates with BH mass 
jFranceschini et alll998llLacv et al.l20o'll : lJarvis fc McLurd 
l2002ft. Also the fun d amental plane of black-hole activity 
jMerloni etafl 120031: iFalcke et all l2004) suggests a weak 
mass dependence as the radio luminosity depends nonlin- 
early on the accretion rate (not scaled to Eddington units). 
However, other researchers have found no evidence for corre- 
lations be tween radio loudness or radio power and black-hole 
mass fsee lHcll2002UWoo fc Urnl2002l e.g.), as selection ef- 
fects severely affect most of the studied samples. Thus, it 
is important to explore whether a mass dependence of the 
radio loudness can modify our findings. 

In the left-hand panel of Fig. we plot the average 
radio loudness against disc luminosity, i.e., the optical lu- 
minosity. We do not find a strong correlation between the 
two as long as we consider all SDSS quasars. Only when we 
exclude X-ray non-detections, a clearer trend is found. But 
also in the first plot the most luminous objects should have 
— on average — larger BH masses than the faint objects. 
Thus, a strong mass effect is not seen for all SDSS quasars 
with radio detections. Furthermore, it is well- known that 
radio -loud objects are more X-ray bright (e.g., lElvis et all 
Il994l) . This supports the idea that the observed dependence 
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of radio loudness on position in the DFLD is not purely due 
to different black hole masses. 

Jet models often assume a linear coupling between the 
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larger black holes should have larger radio luminosities as 
the physical accretion rate is higher for a fixed Eddington 
ratio. In our DFLD we show the radio loudness, i.e. the ratio 
of the radio flux and the flux in the B-band. A constant R 
indicates that the radio flux scales linearly with optical flux 
and probably the accretion rate. However, it may be that the 
radio luminosity scales with a higher power of the accretion 
rate. E.g., for the radio core, the radio luminosity scales 
roughly with L-Rad oc M if the jet power is a constan t 
fraction of the accretion rate telandford fe Konigl Il979l) . 
This scaling is further sup ported by the "fundamental plane 
of accreting black holes" iFalcke et alJl2ofl iMerloni et alJ 
l2003l:lKording et al-llioOrj . We have verified that this effect 
does not change the plot of the radio-loudness in the DFLD. 
The figure looks similar if we plot Ljt M d/L^ 4 (this assumes 
that the disc luminosity is a good tracer of M). Thus, this 
dependence of the radio luminosity on accretion rate cannot 
reproduce the observed diagrams. 

It has been suggested that instead of a correlation be- 
tween R and black-hole mass, there is a hard switch in ra- 
dio loudness for source with BH m asses M > 10 8 M Q (e.g., 
lLaorll200Cl:lMcLure fc Jarvisll2004h . The Eddington limit for 
a 10 s M Q BH is around 10 46 erg s" 1 , but our DFLD still 
shows a gradient in radio loudness above this luminosity. 
This gradient can therefore not be solely due to a single 
switch. Thus, if the BH mass is the dominant effect in the 
DFLD, there must be an additional mass dependence be- 
sides the suggested switch around M > 10 8 Mq. 

As a "worst-case scenario", let us assume that all ob- 
jects have a fixed Eddington ratio. In this case, the lumi- 
nosity range of 2 orders of magnitude is due to 2 orders 
of magnitude difference in the black-hole mass. The highest 
and lowest average values of R found in our DFLD differ 
by a factor 10 2 ' 5 . If the differences in radio loudness were 
purely due to a mass dependence, this would suggest that 
Lr tx A1 2 S . Furthermore, if this effect was to explain the 
entire radio-loudness distribution on the DFLD, the most 
massive black holes would have to be the most X-ray bright 
objects (i.e., lie on the right-hand si de of the diagram ). This 
does not seem to be the case ('e.g.. lWoo et al.ll2005h . Since 
the req uired mass scaling seems to be rather large (compare 
Fig. 2 in iMcLure fc ,Iarvisl2004t) and the needed X-ray to op- 
tical properties are not found, we consider it unlikely that 
the radio loudness distribution in the DFLD arises purely 
due to a mass scaling of the radio loudness. 

4.3 Orientation, beaming, and host galaxy light 

The orientation of a quasar's accretion disc and jet with 
respect to the observer can pl ay an important role fo r the 
observed luminosities (see e.g.. lUrrv fc Padovanill9'95Tl . The 
main effects are relativistic beaming of the emission origi- 
nating from the jet and obscuration of disc and broad-line 
emission by a torus. All our SDSS quasars are of Type 1, so 
broad lines are visible in all objects. It follows that we se- 
lect only sources that are not strongly obscured (see also the 
discussion of the likely fraction of reddened quasars in the 



SDSS bv lRichards et alJl2003l: fHopkins et al.ll2004l) . This is 
especially important as we use soft X-ray data from ROSAT, 
which is very sensitive to absorption. Due to our selection 
of broad-line objects, we ensure that the emission lines are 
not swamped by the featureless continuum of the relativis- 
tic jet (a BL Lac object would not be classified as a broad- 
line SDSS quasar). Hence, we are selecting against strongly 
boosted sources. 

Thus, we can safely assume that the optical emission is 
mostly due to the accretion disc. Radio emission, however, is 
known to be strongly affected by beaming, because it arises 
from a relativistic jet. If radio emission is the only beamed 
component in the SED, any correlation between the strength 
of the beaming and the position of the source in the DFLD 
supports our statement that the jet properties of AGN de- 
pend on both the optical and X-ray properties of the source. 

In case that some or all of the X-ray emission also orig- 
inates from the jet, the X-ray bright sources on the right 
side of the DFLD will preferenti ally be stronger beamed . 
In unified schemes for AGN fe.g.. llJrrv fc PadovanilHoOBl ). 
flat-spectrum radio quasars (FSRQ) are objects whose ra- 
dio and optical continuum emission is enhanced by rela- 
tivistally beamed emission from a jet; hence, part of the 
X-ray emission will als o be enhanced by b eaming. However, 
the "blazar sequence" ijFossati et alJllifflsT ) appears to indi- 
cate that beamed X-ray emission is progressively less impor- 
tant in objects with higher optical (accretion disc) luminos- 
ity (although a number of recent author s have questioned 
the reality of the blazar sequence; e.g., | Ant6n fc Browne! 
l2005HLandt et all20ollNieppola et al.l2006l) . Based on cor- 
relations between X-ray, radio and optical fluxes, a num- 
ber of authors concluded already in the Einstein era that 
there has to be a beamed componen t of X-ray emission in 
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(1995) found a beaming-induced correlation between radio 
and X-ray fluxes bo th for steep-spectrum and flat-spectrum 
quasa rs, while both lGalbiati et alJ d2005l) and lGrandi et alJ 
(2006) report evidence for a beamed contribution in FSRQ, 
but not in other radio-loud quasars. On the other hand, 
lEvans et alJ (120061) found a contribution to the X-ray emis- 
sion of high-power radio galaxies and quasars from both the 
accretion flow (i.e., the disc) and the jet, indicating that the 
X-ray emission is not purely due to a beamed component. 

To consider whether beaming is important in the opti- 
cal, we begin with the observation that the H a equivalent 
width of all types of AGN is essentially constant l|Osterbrockl 
Il989l Fig . 11.6) and independen t of the radio properties 
(see also iBoroson fc Greer] Il992l . § 4.2; this statement is 
not contradicted by the anticorrelation between continuum 
luminosity and equivalent widt h of quasar em ission lines 
known as the the Baldwin effect |Baldwinlll977| . as that ef- 
fect is strongest in high-ionization lines like Ly a and C iv) . 
The constant Balmer-line width implies that beaming does 
not change the optical continuum luminosity significantly. 
Hence, we expect FSRQs to have much higher values of R 
than other AGN purely due to relativistic beaming. At the 
same time, these objects may have also a high non-thermal 
fraction due to beaming of the X-ray emission. Hence, the 
trend for radio-loud sources to be found at higher values of 
the non-thermal fraction may be due to relativistic beam- 
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ing affecting radio and X-ray emission, but not the optical 
emission, in sources with the most powerful relativistic jets 
(the trend of average R increasing to the right in Fig. 
Unfortunately, we do not have spectral information for our 
radio sources, so that we cannot remove the FSRQ from our 
plot and check whether the remaining radio-loud objects still 
show the same trend. 

Nevertheless, it appears unlikely that the observed 
trends are exclusively due to beaming and not connected 
to some intrinsic differences. First, the trend for the av- 
erage R to increase with higher total luminosity cannot 
be caused exclusively by beaming, as there are high- J? ob- 
jects whose total luminosity is more than 50% optical (ther- 
mal disc) emission. Secondly, the radio-quiet and radio-loud 
quasar population c annot be unified using orientation and 
beamin g alone (see iKell ermann et alJll 989l and refere nces 
therein: iRawlines fc Saunderdll99lt iMiller et al.lll993l) . In 
other words, there is an intrinsic difference in addition to 
relativistic beaming that gives some sources a much higher 
R value than most others. We verified that the observed gra- 
dient in average radio loudness on the DFLD is also visible 
if we plot the fraction of radio loud objects (R > 10) on 
the DFLD. Third, there is no correlation between the radio 
and X-ray fluxes for the objects in our SPSS sample , un- 
like in the samples considered bv lGalbiati et ail i2005h and 
iGrandi et alJ (|200^). Hence, even with a beamed contribu- 
tion to the X-rays, the fact that sources with high thermal 
fraction are also more likely to have a high R value also 
has to arise from an intrinsic connection between thermal 
fraction and R other than relativistic beaming. 

In our analysis we use a fixed hydrogen column den- 
sity for all quasars. The used ROSAT fluxes are sensitive 
to changes in the column density. Thus, strong variations in 
the intrinsic column density of the quasars may change the 
position of the sources on the DFLD. However, the majority 
of our sources is at moderate redshifts (z ~ 1) so the ob- 
served photons have been emitted as harder X-rays and are 
less sensitive to intrinsic obscuration. All our SDSS quasars 
have broad lines. According to the "standard unification" of 
AGN, the intrinsic column density is small or negligible in 
these broad line sources (for abs orption see Richard s et alJ 
120031 . e.g., for the column density iGranato et al.lll997l . e.g.'). 
Nevertheless, we cannot rule out completely that absorption 
affects our DFLD. However, for this to become the dominant 
effect, there would have to be a strong anticorrelation be- 
tween radio loudness and hydrogen column density, which 
we consider to be unlikely given that all objects in the sam- 
ple are broad-line objects and selected on optical continuum 
magnitude. 

Finally, neglecting the contribution of host galaxy 
starlight to the SDSS optical quasar magnitudes has the 
tendency to let lower-luminosity and lower-redshift sources 
appear too optically luminous. Furthermore, it artificially 
incre ases the disc fraction and lowers the R value. How- 
ever, IVanden Berk et al.l (120061) have performed a spectral 
decomposition of SDSS quasar spectra and find that only 
20% of the objects considered have a host galaxy fraction of 
more than 50%, and hardly any have a host galaxy fraction 
of more than 90%. Thus, the host galaxy contribution af- 
fects only a small fraction of objects, and only in the lower 
left-hand part of the DFLD. They will not affect the high- 
luminosity sources that create most of the signal in the upper 



right-hand part of the DFLD. Nevertheless, it would clearly 
be desirable to repeat our experiment including host galaxy 
corrections. 

4.4 Size-luminosity evolution of radio sources 

We only have radio data at a single observing frequency, 
1.4 GHz, and with a rather large beam size of 5", corre- 
sponding to physical scales between 17kpc at z = 0.2, our 
low-redshift cutoff, and 42kpc at z — 1.6, where the an- 
gular diameter distance is maximal for our chosen cosmol- 
ogy. This implies that there will be some contribution to 
our "core" luminosities from extended lobe components in 
the FIRST beam. This is not a problem for sources which 
are large enough for the lobes to be clearly separated from 
the core in FIRST imaging, since the SDSS catalog match- 
ing employs a rat her small matching radius (3") for FIRST 
sources. However. iKaiser et al.l il997f) have shown that lobes 
are more luminous when the source is smallest, exacerbating 
the problem for sources where the core cannot be discerned 
from the lobes based on the FIRST imaging. Nevertheless, 
these sources have lobes precisely because they have power- 
ful jets, so that they must have a large core radio flux, too. 
Hence, any contribution from unresolved lobes will have the 
effect of increasing the inferred value of R only in sources 
that have a high radio loudness anyway, but it will not in- 
crease the radio loudness of sources without a jet. In other 
words, lobe contamination makes radio-loud sources even 
more radio-loud, but it does not create false positives in re- 
gions of the DFLD where there are no radio-loud sources. 



5 DISCUSSION: ANALOGIES WITH XRBS 

We have shown in the previous section that quasars are more 
likely to be radio-loud if the X-ray luminosity is large com- 
pared to the thermal disc luminosity and/or the source is 
very luminous. This suggests the radio loudness depends at 
least on 2 parameters: the total accretion rate as well as a 
measure of the relative prominence of the non-thermal X- 
ray component compared to the thermal disc emission. This 
effect is similar to that observed in XRBs, which motivated 
our development of the DFLD. Here we compare the be- 
haviour of AGN with XRBs. 

5.1 Monte-Carlo simulation of an XRB 
population DFLD 

For X-ray binaries, we can study the evolution of single ob- 
jects in the HID during their outbursts (e.g., Fig. lllH . while 
we have to consider AGN as a population. Nevertheless, as 
DFLDs are constructed as an analogue of the HIDs, the over- 
all distribution of a population of XRBs in an HID and the 
distribution of our AGN sample in a DFLD should be simi- 
lar. Indeed, XRBs and AGN populate the right side of their 
diagrams at low luminosities and create a cloud of objects 
in the so ft and IM states fu rther to the left at higher lumi- 
nosities jFender et al.ll2004l) . To compare XRBs and AGN 
in more detail, we simulate a DFLD for an XRB population 
of 100 objects based on our knowledge of the evolution of 
the disc and the power law component in an outburst cy- 
cle of an XRB (a sufficient number of observed HIDs is not 
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Figure 10. DFLD for a simulated sample of XRBs. We simulated 
100 outbursts of an XRB. 



readily accessible to us). We fix th e black-hole mass of the 
simulated XRBs at 10M© (see e.g.. IOrosd|2003l) . 

For each simulated outburst, we first choose the accre- 
tion rate at the hard-to-soft state transition Merit randomly 
from the range 5%-100% Eddington. We use Lr/Ld as mea- 
sure of the radio loudness to retain the same physical mean- 
ing as Lr/Lo for AGN. During its outburst the source moves 
through its canonical states, and we compute luminosities 
for each of the source components using the prescription 
given in Appendix FBI 

Figure E01 shows the simulated DFLD for 100 XRB out- 
bursts. The overall distribution of the simulated XRBs in the 
DFLD is similar to the observed diagram for AGN. Also the 
radio loudness distribution is similar for XRBs and AGN. 
The sources at the right-hand edge of the diagram are the 
most radio-loud objects. In the upper part, where the SDSS 
quasars are located in the AGN DFLD, we find a gradient 
from high average R in top right-hand corner to a low aver- 
age R in the bottom left part, both for the simulated XRBs 
and the observed AGN. For XRBs, the jet is quenched (we 
assumed quenching by a factor 100) in the bottom left cor- 
ner of the upper part of the diagram (Fig. by contrast, 
all sources in the top right-hand corner are extremely radio- 
loud. The gradient in radio loudness is smaller for AGN, and 
we will discuss possible reasons below. At low luminosities, 
both XRBs and AGN are located in the right/power-law 
dominated side of the diagram. 



5.2 Further evidence for common accretion states 
in XRBs and AGN 

The similarity of the DFLDs of XRBs and AGN futher sup- 
ports the idea that AGN have the same accretion states as 
XRBs. In both classes of object, the hard state is character- 
ized by the absence of a strong disc component, a power- 
law component dominating the total luminosity, and the 
presence of a jet. T he jet power is likely to dominate over 
the radiated power jKording et al.ll20od> . In the AGN zoo, 
we classify LLAGN and FR I ra dio galaxies as hard state 
objects, (as already suggested bvlMeierll2 00l|: |Falcke et alJ 
120041: iMaccarone et al.ll2003l : Iilcll2005l: iKording et alJl2006^ . 
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Figure 11. Sketch of the distribution AGN and XRB states 
(scaled to AGN masses) on the DFLD. As discussed in the text, a 
lobe-dominated AGN could undergo a strong change in the core 
properties, moving it from position "1" to position "2", but the 
corresponding change in the lobe properties would be delayed by 
at least the fading time of the lobe, so that it would appear to 
be too radio-loud based on FIRST data. By contrast, AGN at 
position "3" should have no lobes. 



The hard I MS has some disc co ntribution, but the radio jet is 
still active llFender et al]l20o!l) . In the AGN population, we 
suggest that this state corresponds to the radio-loud quasars. 
In the soft IMS and the soft state in XRBs which are even 
more disc-dominated, the jet is not observable in the radio. 
The corresponding AGN class are the radio quiet quasars. 
In Fig. Illl we show a schematic diagram of the distribution 
of the states in our DFLD. 

For XRBs each cell of the DFLD can be associated with 
a given class, as we have assumed that the BH masses of 
XRBs are similar and used a deterministic evolution for all 
sources. Therefore, we find a fast transition in our simulation 
from the radio-loud hard state and hard IMS to the radio- 
quiet states (soft and soft IMS). In our AGN sample, we 
have a large BH mass range and orientation will affect the 
observed luminosities. Thus, the sharp transitions will be 
washed out. Even if a cell was associated, for example, with 
the hard IMS, there would still be a fraction of radio-quiet 
sources in this bin due to different BH masses or source 
peculiarities. 

With long radio observations, it is usually possi- 
ble to find radio emission from all radio-quiet quasars, 
and there are correlations between radio and emission- 
line l uminosity both for radio-quiet and radio-loud objects 
(e.g.. iMiller et al.lll993t iFalcke et al.lll995l : IXu et al.lll99gl : 
iHo fc Pendl200ll) . The difference in radio loudness between 
radio- loud and quiet quasars i s typically around a factor 100 
(e.g., iKellermann et alJ 11989). For XRBs it is often stated 
that the radio jet is quenched (e.g., iFender et alJ 119991) in 
the soft state. However, the upper limits on the radio flux 
on high state objects are not yet low enough to rule out 
radio emission from "quenched" XRBs at a level equivalent 
to that in radio-quiet qua sars (one of the best studied cases 
is XTE J1550-564 where ICorbel et alJ fcOOll) find that the 
radio emission is quenched by at least a factor 50). The anal- 
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ogy between XRBs and AGN therefore suggests that very 
deep radio observations of soft-state XRBs should detect 
radio emission. 

For XRBs we can measure how fast a sources moves 
through the DFLD. For AGN it is not yet known how long 
a source spends in each accretion state. State changes in 
AGN might be triggered b y the same proces s as in XRBs, 
e.g., disc instabilities (cf., iFrank et al.ll20o3) . But already 
different XRBs or different outbursts of the same source 
move through the diagram at different speeds. Furthermore, 
environmental effects can also change the external accretion 
rate dramatically (e.g, during a galaxy merger), thus the 
timescales do not have to scale linearly with black hole mass. 
These timescales will not play an important role for the ra- 
dio luminosity as a function of the position in the DFLD as 
long as the radio emission responds to changes in the accre- 
tion rate on a faster timescale than the timescale of a state 
change. Our radio data is from the FIRST survey, i.e., high- 
angular resolution images at 1.4 GHz. Especially at higher 
redshifts the corresponding emitted frequency is even higher 
and the radio emission is likely to originate predominantly 
from the radio core. The emission region is therefore fairly 
small (parsec-scale) and the radio flux will respond to any 
state changes rather quickly (~ a few years) — similar to 
XRB case (here the response is even faster, ~ hours). How- 
ever, the large-scale radio structure (radio lobes) responds 
much more slowly. It is therefore likely that some sources 
changing from the radio-loud hard intermediate state to 
the radio-quiet soft state still have strong radio lobes, even 
though the core is already in a radio-quiet state. For a 100- 
kpc quasar jet, we see the core change up to 200 ky earlier 
than the near lobe. The importance of this delay depends on 
the sum of the fading time of the lobe and the typical lag be- 
tween the time when we see the core properties change, and 
the time when the lobe receives the information that the core 
properties have changed. On the DFLD in Fig. 1111 we have 
marked three positions with the numbers 1 to 3. In the area 
marked as "1" we expect to observe predominantly sources 
in the hard intermediate state, i.e., sources with an active 
radio core and possibly strong radio lobes. Some sources in 
the area "2" may have just changed their accretion state to 
the soft state. Thus, some sources should have a weak radio 
core but still have strong radio lobes. In comparison, area 
"3" should be populated mainly with sources with no lobes 
and only weak radio cores. 

5.3 Other DFLD trends: timing and jet Lorentz 
factor 

In XRBs, there are correlations between the source's posi- 
tion in the DFLD and other observables than just the radio 
loudness. Looking for their presence in AGN will be a strong 
test of the AGN-XRB unification scenario. 

The timing features of XRBs depend on the accre- 
tion state of the object and therefore on its position in the 
DFLD. Hard-state objects on the right side of the diagram 
show strong fiat-topped band limited noise (see Ivan der Klisl 
2006). Soft-state objects on the lower left part of the DFLD 
(area "3" in Fig. II IB should mainly show weak power- 
law noise. During the IMS strong quasi-periodic oscillations 
(QPOs) are found. Near the transition from the hard to the 
soft IMS, objects often show strong QPOs around 6 Hz (type 



B: IWiinands et al.lll99st iBelloni et al.ll2005l) . With a linear 
mass scaling, this would translate into a period of ~ 200 d 
for objects in area "1" and "2". However, as with the radio 
loudness map, different black-hole masses and source pecu- 
liarities (obscuration, beaming) will "smear out" the sources 
in the DFLD, so there will not be a sharp transition as found 
in a single XRB outburst. 

Also the bulk Lorentz fac tor of XRB jets m ay depend 
on the position in the DFLD dFender et alJl2004l) . It would 
be interesting to see if this is also the case for AGN. To do 
so, one would need a large sample of AGN with a measure of 
their bulk Lorentz factors. However, reliable measurements 
of Lorentz factors are hard to obtain ( Fender 2003). If we ap- 
proach the jet line in the DFLD, the jet in X RBs beco me in - 
stable and one observes rapid ejections jFender et alJl2004f) . 
It is not yet clear whether this behaviour is also found in 
AGN. One caveat is that while hard-state XRBs seem to 
have rather slow jets, it is thought that those in BL Lac 
objects are highly relativistic. According to the orientation 
unification schema of AGN, the parent population of BL 
Lac objects are low-luminosity FR I radio galaxies, which 
are likely to be hard state objects. Thus, it may be that 
the jet velocities are different for stellar and supermassive 
black holes or that the observed slow jets have a fast spine 
that has not yet been observed in XRBs. A detailed study is 
needed to decide whether the speeds of AGN jets are similar 
to those of XRB jets. 



6 CONCLUSIONS 

We have shown that the hardness-intensity diagrams used 
to diagnose the states of XRBs can be generalised to disc- 
fraction/luminosity diagrams (DFLDs), which can be used 
both for AGN and XRBs. The shape of the distribution 
of low-luminosity AGN and SDSS quasars in this diagram 
is similar to that of XRBs. Additionally, supermassive and 
stellar-mass black holes show a similar dependence of ra- 
dio loudness on position in the DFLD. This suggests that 
radio loudness can only be understood in context of a two- 
dimensional diagram, i.e., as function of non-thermal frac- 
tion and luminosity, and is not directly related to a single 
variable like the accretion rate or black-hole mass. The simi- 
larity of the diagrams for AGN and XRBs supports the idea 
that AGN have the same accretion states as XRBs. 

At low luminosities, AGN and XRBs are always in their 
hard state, which is characterized by strong X-ray and radio 
emission compared to their disc emission. At higher lumi- 
nosities, AGN can either be in their radio-loud hard IMS 
(radio-loud quasars) or their soft state or soft IMS (radio- 
quiet quasar). The position of an AGN in the DFLD deter- 
mines the probability with which it is in each of these states. 
Thus, we find the same coupling between accretion flow/disc 
and jet properties in AGN as in X-ray binaries. 

To avoid all evolutionary and selection effects, it will be 
desirable to repeat our investigation with a sample of AGN 
from a single narrow redshift range, but with high dynamic 
range in luminosity and well-determined masses. Ideally, the 
sample would extend down to LLAGN, which can to date 
only be found in the nearby universe, but not in the same 
volume as large numbers of quasars. Further strong tests 
of the AGN-XRB unification and our proposed equivalence 
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between XRB states and members of the AGN zoo will be 
performed by looking for correlations between DFLD posi- 
tion and other properties, such as QPO frequencies. 
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APPENDIX A: SQL QUERY FOR SDSS 
QUASAR SAMPLE 

The query selects all objects from the SDSS spectroscopic 
database that are classified with high confidence as a quasar, 
and, where available, ROSAT and FIRST parameters. 

select . . . 

from specphoto as sp 

left outer join rosat as ro on ro.objid = sp.objid 

left outer join first as fi on sp.objid = fi.objid 

where specclass in (3,4) 

and zconf >= 0.35 

and z between 0.2 and 2.5 



Cyg X-l, e.g., is nearly always near the transition and does 
not show any hysteresis. Thus, we set log(0.1MbackC 2 ) = 
37—0.5 (log(0.1M C ritC 2 )— 37), where all luminosities are mea- 
sured in ergs -1 . For this transition we use a similar prescrip- 
tion for the power-law luminosity as in the first transition: 
Lpl = (1 - C)0.1M bac kC 2 and L D = £ 0.1Af back c 2 , but the 
jet is assumed to be quenched during the whole transition 
and is expected to restart only when the source enters its 
hard state. 

• Declining hard state: same as initial hard state, just 
M cr it is exchanged with Mback 



APPENDIX B: RECIPE FOR SIMULATING AN 
XRB OUTBURST CYCLE 

Here, we describe our assumptions on the evolution of an 
XRB during a single outburst. A source is assumed to be in 
quiescence before the start of the simulation, i.e., we start 
the simulation with the rise of the outburst. 

• Initial hard state: When a source starts its out- 
burst, it is in the hard state with its SED dominated 
by the power law component. Therefore, we set Lpl = 

(^jj^—^j 0.1A/ cr itC 2 , assuming a quadratic scaling of the lu- 
minosity with accretion rate as expect ed for inefficient flow 
or jet mo dels (theoretical pred iction: iNaravan fc Yil fl994L 
empirical: iKording et alJ I2006T) . The jet is active and the 
radio luminosity Lr oc M 1 ' 4 ijBlandford fc Koniell Il979t 
iFalcke fc Biermanr]|l995h . The disc is ususally not directly 
visible in this state, we set its luminosity to 1 % of the PL 
component. 

• IMS: During the transiti on, the bolometric luminosity 
does not change significantly dZhang et alj fl*997l . Thus, we 
set L PL = (1 - x 0.1M crit c 2 and L D = £ x 0.1M crit c 2 , 
where £ is the disc fraction during the transitio n. The radio 
jet is active up to the jet line iFender et al.l20o3) . which may 
vary with the critical accretion rate M cr it . Higher-luminosity 
outbursts have an active jet up to larger disc fractions £. We 
set the jet line at £ = 0.57 + -j^critc j g norma ii sec i 
for th e transition in GX 339-4 observed by iBelloni et alJ 
(2005). The jet power up to the jet line is assumed to be 
constant and equal to the maximum jet power achieved in 
the hard state (i.e., both the radiative and the kinetic lu- 
minosity remain constant during the transition). We omit 
the relativistic ejections which happen right at the transi- 
tion to the radio-quiet phase, as they are very short-lived. 
Their inclusion would yield a higher radio flux right at the 
jet line. The strength and size on the DFLD of this feature 
depends on the timescales of a source moving on the DFLD 
compared to the timescale of a responds in the radio (see 
discussion in sec. 15. 21 . 

• In the soft state, the SED is dominated by the soft- 
spectrum, efficiently radiating accretion disc: Ld = O.lMc 2 . 
A hard power law is often visible up to very high photon en- 
ergies, we assume Lpl = O.OILd- The radio jet is quenched 
by a factor 100 compared to the hard state. 

• IMS: The transition back towards the hard state hap- 
pens at lower luminosity than the hard-to-soft transition. 
The hysteresis is smaller for smaller transition luminosities; 
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